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Abstract 

a: 

Qh. An SO(W) supersymmetric grand unified model proposed earlier leading 

<-j ! to the solar solution involving "just-so" vacuum oscillations is reexamined to 

study its ability to obtain the other possible solar solutions. It is found that 
all four viable solar neutrino oscillation solutions can be achieved in the model 
simply by modification of the right-handed Majorana neutrino mass matrix, 
Mr. Whereas the small mixing and vacuum solutions are easily obtained with 
several texture zeros in Mr, the currently- favored large mixing angle solution 
requires a nearly geometric hierarchical form for Mr that leads by the seesaw 
formula to a light neutrino mass matrix which has two or three texture zeros. 
The form of the matrix which provides the "fine-tuning" necessary to achieve 
the large mixing angle solution can be understood in terms of Froggatt-Nielsen 
diagrams for the Dirac and right-handed Majorana neutrino mass matrices. 
The solution fulfils several leptogenesis requirements which in turn can be 
responsible for the baryon asymmetry in the universe. 



Typeset using REVTeX 



* electronic address: albright@fnal.gov 

t electronic address: smbarr@bartol.udel.edu 



1 



I. INTRODUCTION 



Recent results from the Super-Kamiokande Collaboration [[U involving atmospheric neu- 
trinos have rather convincingly demonstrated the partial disappearance of muon-neutrinos 
and favor the oscillation of muon-neutrinos into tau-neutrinos, rather than into sterile neutri- 
nos at the 99% confidence level. With regard to solar neutrinos, the situation is somewhat 
more ambiguous. On the basis of the recently announced 1258 day sample results from 
Super-Kamiokande 0, together with the flux data from the Chlorine [[| and Gallium ||] 
experiments, the partial disappearance of electron-neutrinos through oscillations into the 
active flavors of muon- or tau-neutrinos is favored over oscillations into purely sterile neu- 
trinos, with the large mixing angle (LMA) solution strongly preferred over the small mixing 
angle (SMA), the LOW, and the quasi- vacuum (QVO) solutions. Several recent analyses 
H based on the smaller 1117 day sample are basically in agreement with this conclusion by 
Super-Kamiokande but assign slightly higher probabilities to the other three solutions than 
does ref. 0. 

Whereas the data at present prefer the LMA solution to the solar neutrino problem, from 
a model building point of view the LMA solution seems by far the most difficult solution to 
obtain ||. Many published models of neutrino masses and mixings either cannot obtain the 
LMA solution, or can only obtain it by fine tuning parameters. It is thus of importance to 
reexamine various approaches to see whether they have sufficient flexibility to accomodate 
the LMA solution in a natural way. 

One approach that is particularly flexible is the so-called "lopsided mass matrix" ap- 
proach. The idea here is that the large atmospheric neutrino mixing angle arises from the 
form of the charged lepton mass matrix. In other words, in this approach is more nat- 
urally thought of as a mixing of \i and r rather than of and v T . On the other hand, the 
solar neutrino mixing can come from the neutrino mass matrix. In this way the atmospheric 
neutrino problem and the solar neutrino problem can be decoupled from each other. This is 
one feature that allows the lopsided mass matrix models to be more flexible in dealing with 
the solar neutrino problem. In this paper we study an especially simple but very predictive 
example of a lopsided mass matrix model to see whether it can accommodate the LMA 
solution in a natural way, that is, without fine-tuning. 

The model we shall discuss was developed in a series of papers 0, @, by the present 
authors, together with K.S. Babu earlier in the collaboration. The model is based on super- 
symmetric 50(10) grand unification. As is well known, 50(10) symmetry typically relates 
the forms of the Dirac mass matrices of the up quarks, down quarks, charged leptons and 
neutrinos (which we denote by U, D, L, and N, respectively) very closely to each other. In 
this model, the lopsidedness of the charged lepton mass matrix, L, and of the down quark 
mass matrix, D, allow an elegant explanation of many of the features of the quark and 
lepton masses and mixings; in particular, the fact that is large whereas V c b is small. An 
interesting point is that in this model the largeness of the atmospheric mixing angle is 
forced upon one by the structure of L, which in turn is tied by 50(10) symmetry to the 
forms of the other Dirac mass matrices. On the other hand, as is again typical of 50(10) 
unification, the Majorana mass matrix of the right-handed neutrinos, M#, is only indirectly 
related to the Dirac matrices, and is therefore much less constrained. This allows various 
possibilities for solar neutrino mixing. 



In the first papers describing this model 0, it was found that the SMA solution is very 
easily obtained if one assumes certain simple forms for Mr, specifically ones which have 
zeros in the 12, 13, 21 and 31 elements. Later it was realized that the QVO solution is also 
easily obtained || by assuming certain other simple forms for Mr. However, it was found 
that the simplest looking forms for Mr, namely those with many texture zeros, cannot give 
the LMA solution || . In light of the recent claim that the LMA solution is strongly favored, 
we re-examine this model to see whether the LMA can be obtained in a natural way. In 
fact, we look at all four solar solutions. 

In Sect. II we specify the conditions for each of the four solar solutions. The Dirac mass 
matrices and parameters obtained earlier for the SO (10) model in question are presented in 
Sect. Ill, where we also numerically determine the structures of the right-handed Majorana 
matrix needed to reproduce all four solutions. A survey of these numerical results in Sect. IV 
reveals that Mr for the LMA solution, in particular, has a remarkably simple texture which 
can be easily related to the Dirac neutrino matrix. For this case, the seesaw mechanism then 
leads to a light neutrino mass matrix which has two or three texture zeros. The implications 
of this solution for leptogenesis are briefly discussed. 



II. PREFERRED REGIONS IN THE NEUTRINO MIXING PLANE 

Here we summarize the preferred points in the neutrino mixing plane for the atmospheric 
neutrino and the four viable solar neutrino oscillation solutions. We use this information to 
reconstruct the Maki-Nakagawa-Sakata (MNS) |K| neutrino mixing matrix for each of the 
four solutions. 

For the atmospheric neutrino oscillations, the best fit values obtained are fl[] 

AmL = 3.2 x 1(T 3 eV 2 , 

(1) 

sin 2 2^ 2 3 = 1.000, 

in terms of Am| = m 2 — m 2 with sin 2 a4m = ^\Ufj,3\ 2 \U T s\ 2 expressed in terms of the MNS 
leptonic mixing matrix elements. Note that to a high degree, the atmospheric neutrino 
mixing is observed to be maximal. The best fit values for the four solar neutrino solutions 
according to an earlier analysis by Gonzalez- Garcia [llj are 



SMA : Am 2 21 = 5.0 x 10~ 6 eV 2 , 

sin 2 20 i2 = 0.0024, tan 2 12 = 0.0006, 
LMA : Am 2 21 = 3.2 x 10" 5 eV 2 , 

sin 2 20 12 = 0.75, tan 2 12 = 0.33, 
LOW : Am 2 x = 1.0 x 10~ 7 eV 2 , ^ 

sin 2 20 12 = 0.96, tan 2 12 = 0.67, 
QVO : Am 2 21 = 8.6 x lO" 10 eV 2 , 

sin 2 20i2 = 0.96, tan 2 12 = 1.5. 

In general the MNS mixing matrix, analogous to the CKM quark mixing matrix, can be 
written as 
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C12C13 



Umxs = | -S12C23 - Ci2S 2 3Si 3 e 

S12S23 - Ci 2 C 2 3Si 3 e lS 



id 



S12C13 

C12C23 — Si2S2 3 Si3e l 
-C12S23 - Sl2C 2 3Sl 3 e 



id 




(3) 



in terms of c\i = cos #12, S12 = sin# 12 , etc. With the oscillation parameters relevant to the 
scenarios indicated above, to a very good approximation 6 13 = 0° and #23 = 45° whereby 
Eq. ([3D becomes essentially real and of the form 



U 



MNS 



( C12 _ S12 
-S12/V2 c 12 /V2 1/V2 



(4) 



where the light neutrino mass eigenstates are given in terms of the flavor states by 



^3 



73 K + ^ 



v 2 = v e sin ^12 + -^(^n - v T ) cos 012, 
v x = v e cos9 12 - - v T ) sin 012- 



(5) 



For the SMA solution, ^12 = 1.4°, while the three large mixing solar solutions differ from 
maximal in that the angle is approximately 30° for the LMA, 39° for the LOW, and 51° for 
the QVO solutions. Numerically we find for each case 



U 



(SMA) 
MNS 



u 



(LMA) 
MNS 



U 



(LOW) 
MNS 



U 



(QVO) 
MNS 



( 0.9997 
-0.0170 

V 0.0170 

/ 0.866 
-0.354 

V 0.354 

/ 0.774 
-0.448 

V 0.448 

0.633 
-0.548 
0.548 



0.0241 
0.7069 0.7071 
-0.7069 0.7071 

0.500 

0.612 0.70 

-0.612 0.70 

0.633 

0.547 0.70 

-0.547 0.70 

0.775 

0.447 0.707 

-0.447 0.707 



(6) 



III. MODEL MASS MATRICES AND NUMERICAL DETERMINATIONS OF M R 

The model we are studying here is an S'O(IO) grand unified model. For details of its field 
content, the flavor symmetry U(l) x Z2 X Z2, couplings, and so forth, the reader is referred 
to the series of papers in which the model was developed R] , [H] ■ Here we will only mention 
a few of the features of the model important for the present considerations. 

This model arose from an attempt to construct a realistic S'O(IO) model with the simplest 
possible, or "minimal," Higgs content. This attempt led very naturally to the following 
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structures at the GUT scale for the Dirac mass matrices of the up quarks, down quarks, 
neutrinos, and charged leptons, labeled U, D, N, and L, respectively: 
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N 



( V 
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5 S'e 1 * 
D = | 5 a + e/3 | M D , 
tf'e** -e/3 1 

(7) 

/ 5 S'e^ 

5 -e | M D . 
\8'e^ a + e 1 





~ 113 GeV, 




~ 1 GeV, 


a 


= 1.78, 


€ 


= 0.145, 


8 


= 0.0086, 


8' 


= 0.0079, 


<t> 


= 54°, 


V 


= 8 x 10- 



A crucial point is that the four Dirac matrices are closely related to each other by the group 
theory of 5O(10) and that their forms are definitely fixed in terms of a few parameters. 
As a result the model is very predictive, and in fact gives excellent agreement with all the 
known facts about the CKM mixings, the quark masses, and the charged lepton masses. By 
fitting these data, taking into account the renormalization effects from the GUT scale to 
low energies, the following values of the parameters were obtained: 



(8) 



A critical feature of the model is that the parameter a is of order unity, and appears in 
an asymmetrical or "lopsided" way in L and D. This fact plays many roles in the model 
and is indeed the key to its economy and success in fitting the data. It explains (a) why 
m c /m t <C m s /mfo, since m c /m t ~ e 2 , while m s /m{, ~ ea; (b) why the Georgi-Jarlskog 
relation m s /m^ = ^m^/m T holds, since without the a term a factor of | rather than | 
would result; and (c) why V c b U^. The reason for the last is that a appears in the 
32 element of L, where it causes a large mixing of left-handed muon and tau leptons, i.e. 
large U^, whereas it appears in the 23 element of D, where it causes a large mixing of 
right-handed quarks, which is not relevant to V^. The mixing V c b is instead controlled by 
the 32 element of D, which is the small parameter e/3. The fact that a appears transposed 
between D and L has to do with the SU (5) structure of the fields involved. 

For present purposes the most important fact is that the largeness of the atmospheric 
neutrino mixing angle comes from the parameter a in the charged lepton mass matrix L. 
The contribution of the neutrino mass matrix to this mixing is formally of order e, as can 
be seen from the form of N, and is therefore numerically small for generic choices of Mr. 
On the other hand, one sees that the solar neutrino mixing angle receives only a small 
contribution from the charged lepton sector, since the 12 and 21 elements of L are small. 
Therefore, whether the solar angle is large or small is controlled by the neutrino mass matrix 
M v = -N T M^N, or in other words by Mr, since iV is fixed. The form of Mr is rather 
independent of the forms of the Dirac matrices given in Eq. (0) because it comes from 
completely different operators. That is why in this model — and indeed in the general 
framework of "lopsided mass matrix models" in which the atmospheric angle arises from 
large lopsided entries in L — there is great flexibility in how the solar neutrino problem is 
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solved. Different solar oscillation solutions can be obtained by changing the form of Mr 
without affecting in any way the fits to the CKM parameters, the masses of the quarks and 
charged leptons, or the fact that the atmospheric neutrino angle is large. 

In our first papers where this model was discussed, forms of Mr were assumed in which 
the SMA solar solution was naturally obtained. Indeed, one sees immediately that if Mr has 
vanishing 12, 21, 13, and 31 elements, M v does not contribute to the solar neutrino angle, 
which then comes entirely from L and is therefore small. 

The QVO solution can also be very easily obtained. In || the following simple form of 
Mr was constructed: 



Mr = Ae 3 A R . (9) 




With this form the seesaw formula [13] gives the light neutrino mass matrix to be 



M u = N T M^N 



I 

e 2 e \M%/A R . (10) 

\-r]/(Ae 2 ) e 




With Ar = 2.4 x 10 14 GeV and A = 0.05, a fairly reasonable fit to the quasi- vacuum solution 
then emerged with 

m 3 = 54.3 meV, m 2 = 59.6 [ieV, mi = 56.5 fieV, 
M 3 = 2.4 x 10 14 GeV, M 2 = M X = 3.66 x 10 10 GeV, 
U e2 = 0.733, C/ e3 = 0.047, = -0.818, 8' CP = -0.2°, 
Am§ 2 = 3.0 x 10~ 3 eV 2 , sin 2 2B atm = 0.89, 
Am 2 x = 3.6 x 10- 10 eV 2 , sin 2 26 so i ar = 0.99, 
Am 2 x = 3.0 x 10~ 3 eV 2 , sin 2 29 reac = 0.009. 

We now wish to search for right-handed Majorana mass matrix textures which fit more 
accurately each of the four solar neutrino solutions. We first note that the MNS mixing 
matrix corresponds to the product of two unitary transformations, 

Umns = U ] L U U , (12) 

where Ul diagonalizes the Hermitian lepton matrix L^L, and U u diagonalizes the light neu- 
trino mass matrix which we assume to be real and symmetric for simplicity: 

L diag tL diag = Jj\p L U u M diag = jjT M jj^ ^3) 

It is easy to see that, given a specific pattern of neutrino masses and mixings, one can invert 
to find a form of Mr that will give that pattern. To be given a pattern of neutrino masses 
and mixings means that one is given the MNS mixing matrix Umns and the neutrino mass 
eigenvalues mi, m 2 , and m 3 . On the other hand, the model itself specifies the charged 
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lepton matrix, L, and the neutrino Dirac mass matrix, N\ cf. Eq. (|7|). Thus Mr can be 
inferred as follows. First, Ul can be directly obtained from diagonalization of L^L. Then 
Ul together with the given Umns determine U v through Eq. ([12]). Although L and hence 
Ul are complex, we can obtain a real U v by making use of the freedom to perform a phase 
rotation on Umns, so that 

U V = U L - diag(l, 1, e"*) • U MNS . (14) 

Then, defining 

ag = diag(m h -m 2 , m 3 ), (15) 

with hierarchical masses chosen which are related to the Am? 's, one can use this and the 
matrix U v already found to determine M v by using the second of Eq. ([13]) . Finally, one can 
use the N known from the model and M v to find Mr by inverting the see-saw formula 

M R = NM~ 1 N T . (16) 

We present the numerical results for each of the four solar solutions as follows: 

/ 0.156 x 1(T 7 -0.190 x 1(T 4 0.116 xl0" 3 \ 
m (sma) = _ 19Q x 1Q -4 o 0105 -0.123 x 5.2 x 10 14 GeV, 

V 0.116 x 10~ 3 -0.123 1.000 / 

with M 1 = 3.7 x 10 6 , M 2 = 2.3 x 10 12 , M 3 = 5.3 x 10 14 GeV; 

/ 8.30 x 10~ 10 -0.511 x lO" 5 2.13xl0- 5 \ 
m (lma) = _ 511 x 10 -s Q 0244 -0.155 x 3.0 x 10 14 GeV, 

\ 2.13 x 10~ 5 -0.155 1.000 / 

with Mi = 4.2 x 10 6 , M 2 = 6.7 x 10 10 , M 3 = 3.1 x 10 14 GeV; 



M (LOW) 



I 5.15 x 10~ 10 -1.43 x lO" 5 5.46 x 10" 5 
= -1.43 x 10" 5 0.0292 -0.176 I x 5.8 x 10' 1 GeY 

V 5.46 x 10" 5 -0.176 1.000 

with Mi = 6.0 x 10 6 , M 2 = 9.7 x 10 11 , M 3 = 6.0 x 10 14 GeV; 
-6.98 x 10~ 10 -1.33 x 10~ 5 4.75 x 10" 5 



(17) 



r (QVO) _ | _ x 33 x 10 _ 5 0.0481 -0.222 1 - * ° - inl5 

4.75 x 10~ 5 -0.222 1.000 

with Mi = 8.8 x 10 6 , M 2 = 3.8 x 10 12 , M 3 = 3.3 x 10 15 GeV. 

Strictly speaking, the above results were obtained at the GUT scale, but with the moderate 
value of tan (3 ~ 5 preferred by the model || and for the hierarchical and sign choices given 
in Eq. (|l^) above, the evolutions in masses and mixings from the GUT scale to the low 
scales are extremely small and can be neglected [|T4|]. 

That one can find forms for Mr that reproduce the various solar neutrino solutions is in 
itself not very significant, for as we have just seen, this is guaranteed as long as the relevant 
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matrices are invertible. The significant question is whether the matrix M R that gives a 
certain solar solution is obtainable in the model under discussion in a simple way without 
fine-tuning. The forms for M^ 5M ^ and M R ® V °^ given in Eq. (|iT|) are complicated-looking. 
However, these are the forms that reproduce the present best-fit SMA and QVO solutions 
according to JTTJ . One already knows from our previous work, as has already been mentioned, 
that much simpler forms for M R , having several texture zeros, give perfectly satisfactory 
SMA and QVO solutions; moreover, those simpler forms are obtainable straightforwardly 
without fine-tuning. But that same earlier work shows that forms for M R having several 
texture zeros do not yield a satisfactory LMA solution in this model. The question is then 
whether the form M R LMA ^ given in Eq. (|17D, or something sufficiently close to it, can be 
obtained simply and naturally in the model. To this question we now turn. 



IV. SIMPLE ANALYTIC FORM FOR M R INVOLVING THE LMA SOLUTION 

At first glance the form of M R LMA ^ in Eq. ([T7| ) looks very complicated. However, it has 
some significant features that suggest that it may be obtainable in a simple way. First of 
all, one sees that (M R ) 2 3 = (Mr) 32 — ~ e and (M R ) 2 2 — e 2 , where e is the parameter that 
appears in the Dirac matrix, N, of Eq. (0). To a good approximation we can therefore 
introduce the analytic form 

M R LMA ^=\-h (ij r* -r|A w . (18) 




written in terms of parameters appearing in the Dirac neutrino matrix, where e = 0.145 and 
T) = 0.8 x 10" 5 as before, cf. Eq. (|), and A R = 2.5 x 10 14 GeV. It will turn out that the 
new parameters a, b and c are all of order unity in order to obtain the LMA solar solution. 
Making use of the seesaw formula, we then find 

/ e/{a-b) \ 

M {LMA) _ e/(fl _ 6) _ e 2 (c 2 _ tfy( a _ b y _ k/ /( a _ 6 ) My/ A R . (19) 

V -be/(a-b) 1 / 

It is interesting, and we shall see, relevant to leptogenesis that this form has some texture 



zeros. These texture zeros follow directly from the form of the 23 block of Eq. ( 18f) . That 



this 23 block has rank 1 immediately suggests that it can arise from diagrams of the Froggatt- 



Nielson type |15| . Moreover, the fact that the same parameter e appears in both M R and N 
suggests the possibility that the hierarchies in the 23 blocks of both matrices may have the 
same origin. These suggestions can be realized as we now show. 

In Fig. 1 we repeat for clarity the diagrams in our model which contributed to the Dirac 
matrices in the 2-3 sector. The dominant 33 elements arise from the 10# Higgs electroweak 
doublet contributions. For the 23 and/or 32 elements, higher-order contributions arise from 
electroweak doublets in both the 10# and 16# 5*0(10) representations, with additional 
singlet Higgs VEV's and a 45# Higgs GUT scale VEV pointing in the B — L direction. Due 
to the SU(5) structure of the Higgs fields, the diagram appearing in Fig. 1(c) contributes 



S 



only to the D 2 s and L 32 elements of the down quark and charged lepton mass matrices. Note 
that the internal superheavy fermions appearing in 16, 16, 10i and IO2 are integrated out. 

In Fig. 2 we show the lower-order diagrams which can contribute to the 2-3 sector of 
the right-handed Majorana mass matrix. Here a singlet Higgs GUT scale VEV, Vm, couples 
two superheavy conjugate singlet fermions thus inducing a breaking of lepton number. The 
VEV's in the 16^'s also appear at the GUT scale. The 16h — 1'h P a h appearing in insertion 
"A" of Fig. 2 serves to lower the heaviest right-handed Majorana neutrino mass down to 
A R = 2.5 x 10 14 GeV from the GUT scale value of 2 x 10 16 GeV. By making use of the 
techniques spelled out in detail in ||, one can readily show that the 23 elements of N 
and Mr are scaled by the same factor e relative to their 33 elements. The factor enters 
antisymmetrically in TV for the 23 and 32 elements due to the B — L nature of the 45# VEV 
and the presence of both left-handed neutrino and conjugate neutrino states, while it appears 
symmetrically in Mr since both states involve conjugate neutrinos. In the Majorana case, 
both superheavy singlet and 45 fermions must be integrated out. We have checked that the 
these diagrams can be achieved as indicated with proper assignment of the U(l) x Z2 x Z2 
flavor quantum numbers for the new heavy fermion fields introduced. 

We now turn to the small entries of the first row and column of M R LMA ^ in Eq. (|I8D with 
a, b, and c numbers of order unity. The fact that the whole matrix manifests a geometrical 
hierarchy involving the same small parameters e and 77 that appear in N reinforces the idea 
that Mr may be simply obtained by Froggatt-Nielsen-type diagrams involving some of the 
same VEVs that generate N. If it were the case that a = b = c exactly, then the whole 
matrix would have rank 1, and thus all its elements could be obtained from a single Yukawa 
vertex (lgl^VM, in the same way that we illustrated for the 23 block. However, that would, 
of course, be unrealistic in that two neutrinos would then be massless. However, it is not 
necessary that the matrix be of rank 1 in order that it arise from simple Froggatt-Nielson 
diagrams. Thus we have the possibility that a, b, and c are not all equal. 

For an especially interesting numerical example, suppose that 



This has a simple interpretation in that all elements of the Mr matrix receive contributions 
from the Yukawa vertex involving Vm, while only the 13 and 31 elements receive contributions 
from a second AL = 2 violating Yukawa vertex involving V' M . This can be realized with the 
proper choice of flavor indices for V' M . By the see-saw formula, one then has 



a = 1, b = c = 2, A = 2.5 x 10 14 GeV. 



(20) 




(21) 



with three texture zeros from which we obtain 
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?»3 = 57.4 meV, rri2 = 9.83 meV, m\ = 5.61 meV, 



M 3 = 


2.5 x 10 14 


GeV, M 2 


= Mi = 2.8 x 10 8 GeV, 


u e2 = 


0.572, U e3 


= -0.014 


= 0.733, 5 CP = 0° 


Am 2 2 


= 3.2 x 10 


" 3 eV 2 , 


sm 2 26 atm = 0.994, 


Am 2 2l 


= 6.5 x 10 


" 5 eV 2 , 


sin 2 2fl soiar = 0.88, 


Am 2 31 


= 3.2 x 10 


" 3 eV 2 , 


sin 2 29 reac = 0.0008. 



(22) 



These results fit both the atmospheric and the LMA solar mixing solutions extremely well 
and can be considered a success for the model. In fact, the best fit point for the LMA solar 
mixing solution as given by the Super-Kamiokande Collaboration in their latest analysis 
of 1258 days of data taking @ is (sin 2 2# so ; = 0.87, Am 2 21 = 7 x 10~ 5 eV 2 ). We find the 
whole newly-allowed LMA region can be covered with a, b and c varying by factors of 0(1) 
from the values given in Eq. (|20|) . It is noteworthy that the solar neutrino mixing is near 
maximal, but not actually maximal as that is presently excluded experimentally by the 
SuperKamiokande results at more than the 95% confidence level. 

How fine-tuned is the form of Mr that we have been discussing? One feature that at 
least appears fine-tuned is the fact that the 23 and 32 entries in Eq. (|T8| ) are not only of 
order e but actually equal to — e exactly. This one has no right to expect from the mere fact 
that the same VEVs come into the diagrams for N and Mr, since as can be seen from Figs. 1 
and 2 different Yukawa couplings are involved in the 23 entries of the two matrices. One can 
test how fine-tuned the form in Eq. ( |TS| ) is by replacing the 23 and 32 elements by -de and 
the 22 element by d 2 e 2 . (The fact that the same d enters is due to the factorized structure 
of the diagrams in Fig. 2, and is therefore not a fine-tuning.) One naturally expects that d 
is of order unity, but how close must it be to 1 to give a realistic LMA solution? It turns out 
that the most severe constraint on the value of d comes from the limit on U&. To satisfy the 



condition that \U&\ < 0.15 |T6[, one requires that 0.85 < d < 1.15. Thus, the LMA solution 
does not require an unnatural fine-tuning of parameters. 

Finally we note that the upper bound on the lightest heavy Majorana neutrino mass 
Mi should be less than or of order 10 9 GeV to prevent overproduction of gravitinos from 
overdosing the universe after inflation |T7J . This bound is satisfied for all four solar solutions 



as determined in Eq. ( pTTD and, in particular, for the model illustrated above. A second 
condition for leptogenesis is that the 13 and 31 elements of Mr be suppressed by a factor of 
at least 10 3 relative to the 33 element to inhibit mixing of the heaviest right-handed neutrino 
with the lightest one in order to prevent its rapid decay washing out the lepton asymmetry 
generated. This is satisfied in our model. 



V. SUMMARY 

We have investigated how an SO(10) SUSY GUT model proposed earlier can be modified 
in order to obtain solar neutrino solutions other than the vacuum solution. The study 
revealed that only the right-handed Majorana neutrino mass matrix needed to be modified, 
with the Dirac matrices for the neutrinos and charged leptons (as well as for the quarks) left 
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unchanged. In short, in this model the maximal atmospheric neutrino mixing is controlled 
primarily by the structure of the charged lepton mass matrix, while the type of solar neutrino 
solution is largely determined by the form of the right-handed Majorana mass matrix. 

Of particular interest was the finding that the large mixing angle solar solution is readily 
obtained with a nearly geometrical hierarchy in M R , where the 2-3 subsector has a close re- 
lationship with that for the Dirac neutrino matrix, as seen by study of the Froggatt-Nielsen 
diagrams. It is precisely this structure which provides the "fine-tuning" necessary to achieve 
the LMA solar solution. 
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FIG. 1. Diagrams that generate the 33, 23 and 32 elements in the quark and lepton Dirac mass 
matrices shown in Eqs. (0). (a) The 33 elements denoted "1". (b) Antisymmetric contributions 
denoted "e," to the 23 and 32 elements where the VEV of the 45# appears in the B — L direction, 
(c) Asymmetric contributions to the 23 and 32 elements denoted "a" appearing in the down quark 
and charged lepton mass matrices arise from this diagram. They do not appear in the up quark 
and Dirac neutrino mass matrices due to the SU (5) structure of the fields explicitly indicated in 
the diagram. 
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FIG. 2. (a), (b) and (c), respectively, show the diagrams leading to (Mr)^, (Mr)^ = (Mr)23 
and (Mp>) 2 2- Note that these diagrams all come from the same vertex (lgl^V^ and so lead to 
an exact factorized or geometrical form where (Mn)^(M r)22 = (-^jj)32(-^,r)23- The insertions 
denoted "A" and "B" are defined in (d) and (e). The ratio B/A is proportional to (45#)/(l" h) 
and so is of order — e = N23/N33, as can be seen by inspecting Fig. 1(a) and (b). 
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